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The Staggered 77' with Smeared Operators 
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We present a refined calculation of the 7/ mass using staggered fermions and Wuppertal smeared operators 
to suppress excited state contributions. We use quenched and dynamical configurations of size 16 3 x 32, with 
Nf = 0, Nf = 2 and Nf = 4, and compare our results with the expected theoretical forms from quenched, 
partially quenched, and unquenched chiral perturbation theory. 



1. INTRODUCTION 

The pseudoscalar spectrum of QCD consists of 
an octet of mesons which are approximate Gold- 
stone bosons of spontaneously broken SU (3) axial 
flavor symmetry, plus an anomalously heavy fla- 
vor singlet meson, the 77'. The heaviness of the 77' 
is attributed to the effects of topology yj . In an 
SU(3) symmetric world, m v > would obey 

22,2 / -1 \ 

m v , = 777.Q -f" 777g (1) 

where m\ is the average mass squared of the octet 
mesons which vanishes in the chiral limit, while 
777,0 i s the topological contribution which does not 
vanish in the chiral limit. Neglecting 77-77' mix- 
ing, one plugs in the known meson masses to 
obtain the "experimental" value mo{Nf = 3) = 
860MeV. 



2. LATTICE CALCULATION OF 

Extraction of m^i from first principles has re- 
ceived much attention . The focus in all these 
studies has been to calculate the ratio R(t) de- 
fined below 



R(t) 
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2— loop 
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where rf if) is the operator that creates or de- 
stroys an 77' meson (in terms of quark fields, for 
staggered fermions, this becomes Qj5®IQ). Two 
point correlation function of this operator yields 
both a disconnected diagram (referred to as 2- 
loop in cqn 0) and a connected diagram (1-loop). 



On dynamical configurations, when m va i — 
rridyn, R(t) takes the following asymptotic form 
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(3) 



where N v is the number of valence fermions, Nf is 
the number of dynamical fermions, B is a constant 
and Am — m^ — trig,. 

For the quenched configurations, infinite itera- 
tion of the basic double pole vertex does not exist 
and it can be shown that the ratio is a linear func- 
tion of time 



R(t) = const. H - t 

2ms 



(4) 



This style of calculation was employed by the 
authors of || who obtained a result of rno(Nf — 
3) = 751(39)MeV using quenched configurations 
and Wilson fermions. We used staggered fermions 
and both dynamical and quenched configura- 
tions and reported a value of mo(Nf = 3) = 
730(250)MeV extracted from dynamical config- 
urations in 

3. SIMULATION DETAILS 

The parameters of the ensemble used in the 
simulation are shown in table jy. For all of the 
configurations listed in table ^ the inverse lat- 
tice spacing is about 2 GeV as obtained from 
m p . m va i for the quenched configurations has 
been chosen 10% higher than that corresponding 
to dynamical {md yn — 0.01, Nf — 2) so that m% 
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is same for both. Propagators were computed 
using conjugate gradient on the 128 node OSC 
Cray T3D machine. For details concerning per- 
formance, the type of source, the method adopted 
for calculating the disconnected propagator etc., 
the reader is referred to 101 . 
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The Statistical Ensemble 
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3.1. Smearing 

The disconnected data is noisy and clean sig- 
nals exist only for the first few time slices. Smear- 
ing helps in reducing excited state contributions 
from the first few time slices, thus enabling a re- 
liable extraction of Am. We use the gauge in- 
variant Wuppertal smearing technique. For the 
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Figure 1. Correlators from Smeared Operators 

connected contraction we computed two propa- 
gators, one with a smeared source and point- like 



sink and the other with a point-like source and 
smeared sink (see Fig. |l|). Correspondingly, for 
the disconnected loops, the sink end was smeared 
in the same way(SS). 
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Figure 2. Ratio from SS and LL correlators at 
m va i = m dyn = 0.01 

Fig. H compares the ratio plot with and with- 
out smearing on dynamical configurations. In the 
initial few time slices both the data are different 
but they begin to coincide after a few time slices 
as they should. 

4. RESULTS 

Fig. H shows the form of the ratio on all the 
configurations listed in Table [j]. For Nf — 2 and 
Nf — 4, the points shown are for m va i = mdyn, 
while for Nf = we plot m va i — 0.011. It is 
gratifying to see that this observable clearly dis- 
tinguishes the number of flavors. 
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Figure 3. Nf dependence of ratio 

One extracts from a linear fit to the 
quenched data. Fig. |J shows m\ versus m va i ob- 
tained from both local and the smeared data, fit 
linearly and extrapolated to the chiral limit. As is 
expected, m\ does not vanish in the chiral limit. 
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Figure 4. Chiral extrapolation of m§ 

The Nf — 2 data shown in Fig. || is fit to the 
form of eqn ^| From the fit one extracts Am and 
hence m^i for all the values of mdyn shown. It 
can be seen in Fig. || that m v ' does not vanish in 
the chiral limit. 



set to unity at lowest order, the data prefer val- 
ues 20-30% larger, an amount which could eas- 
ily be accommodated by higher order chiral and 
0(1/N C ) corrections. 
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Figure 6. Z'/Z versus rudyr, 




Figure 5. Chiral extrapolation of mi, 

When m va i ^ md y m the ratio R(t) takes the 
form 

R{t)= -f [At - Bcxp(-Amt) + C] (5) 

Our data are not precise enough to allow a four 
parameter fit, but using lowest order PQ\PT 
(and neglecting the small momentum dependent 
self-interaction a) we can express A, B and C 
in terms of one unknown parameter, m\. For 
m va i > m dy n (N f = 2, m dyn = 0.01) data we 
get reasonable fits and find a partially quenched 
result of m n i (N f = 3) = 876 ± 16 MeV, remark- 
ably consistent with the fully quenched and fully 
dynamical data. For m va i < mdyn the x 2 is n °t 
reasonable. On the other hand, one need not ven- 
ture far from lowest order PQxPT to fit the data. 
For example, from two parameter fits to the dy- 
namical data, one obtains Z'/Z, the ratio of the 
residues for creating 77 and rj' . While this ratio is 



5. CONCLUSIONS 

The value, we obtain for m n i in the chiral limit 
are summarized in the table below. Smeared 
operators have produced a lower value for m n i 
than that obtained from local operators. Within 



Table 2 

(Nf = 3) in the chiral limit 



Quenched 


Dynamical (m dyn = m val ) 


(MeV) 


(MeV) 


LL 1156(95) 


974(133) 


SS 891(101) 


780(187) 



the quoted statistical errors, the values obtained 
above are consistent with experiment. 
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